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1 Introduction 

There is no need to emphasize the importance of finding rehable non-perturbative 
calculational techniques that can be apphed to strongly interacting theories 
such as technicolor. Without such methods, strongly interacting theories 
can never be tested since no quantitative predictions that can be compared 
directly with experiment can be made. 

The analytic continuation by duality (ACD) technique, which was first 
developed in Ref. []I|, has been proposed in Ref. as a potentially rehable 
way to compute the oblique correction parameter 5* for technicolor theories. 
The advantage of the ACD technique was that it could be applied to both 
QCD-like and walking technicolor theories whereas the dispersion relation 
technique used by Peskin and one of us in Ref. 0] could only be applied to the 
former. Furthermore, the ACD estimate of S for walking technicolor implied 
that walking dynamics could render S negative, making it compatible with 
the current experimental limit. This was in contrast to the result of Harada 
and Yoshida who used the Bethe-Salpeter equation approach to conclude 
that S was positive even for walking theories. 

However, the robustness of the ACD technique has been questioned in 
the literature |^. In this letter, we further investigate the reliability of the 
ACD technique. In section 2, we review the definition of the S parameter 
and explain the ACD technique. In section 3, we apply the ACD technique 
to the perturbative spectral function that represents the contribution of a 
degenerate doublet of heavy fermions to see if the well known result I/Gtt 
could be reproduced. In section 4, we apply the ACD technique to a model 
function which has an imaginary part more representative of the actual QCD 
spectrum. Section 5 concludes. 

2 The ACD technique 

The ACD technique consists of using the analyticity of the vacuum polar- 
ization function to convert the dispersion integral for S, which is an integral 
along the real s axis, into an integral around a large circle in the complex 
s plane. The value of the integral on the circle is then estimated by the 
Operator Product Expansion (OPE) [^. 

We denote by IIxy the g^" part of the vacuum polarization tensor between 



,2 



two currents Jx and Jy'- 

The 5" parameter is defined in Ref. [Q as 

s = -An [nw(0) - num 

where the subscripts V and A denote the neutral isospin-1 vector and axial- 
vector currents, respectively, and the prime denotes a derivative with respect 
to s. 

For latter convenience, we define the functions J-'y{s) and J^a{s) as 

nvv(s) = sj'v(s), 

nAA(s) = IIaa{0) + sJ^a{s), 

and the following shorthand notation for the difference between J-'y{s) and 
^a(s): 

Then, 

S = -47r[^v(0) - ^a(O)] = -47r^(0). 

Note that our notation is slightly different from either Ref. or Ref. 0] so 
care is necessary when comparing formulae. 

The function J^{s) = J-'y{s) — J^a{s) is analytic in the entire complex s 
plane except for a branch cut along the positive real s axis starting from the 
lowest particle threshold contributing to J-'is). Therefore, applying Cauchy's 
theorem to the contour C shown in Fig. |l|, we find 



S = -47rJ^(0) 
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If the radius of the contour R is taken to infinity, the integral around the 
circle a.t \s\ = R vanishes and we obtain the dispersion relation 

S = -4 f dsl^!^^, 
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Figure 1: The contour C which avoids the branch cut along the real s-axis. 



which was used in Ref. ^ to calculate S. However, the dispersion relation 
approach requires the knowledge of InijF(s) along the real s axis which is 
only available for QCD-like theories. 

The basic idea of the ACD technique, on the other hand, is to approximate 
the kernel 1/s by a polynomial 



N 



-^PN{s) = Y.^n{N)s'', SE[So,R] 



n=0 



and use it to make the integral along the real s axis vanish instead. Applying 
Cauchy's theorem to the product pn{s)J-'{s) over the same contour C yields 







-Air 



1 /•« 1 / 

— / dsp]y(s)lm.J^{s) + -— f dsp]\f(s)J^(s) 

IT Jso 2m J\s\=R 



(2.2) 



Subtracting Eq. (p.2|) from Eq. ( |2.1|) , we obtain 

S = Sn -\- Afit, 
where 



S]y = f ds 

I J\s\=R 
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fit 
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ds 



--PN{S} 



InijF(s) 



For sufficiently large N, we can expect Agt to be negligibly small. In fact, it 
converges to zero in the limit N ^ oo (though how quickly the convergence 
occurs depends on the interval [so,-R]). We can therefore neglect it and 
approximate S with S^ which is an integral around the circle \s\ = R only. 
We call Afit the fit error. 

In fitting Pn{s) to 1/s, the s-coefficients a„(A^)'s are determined so that 
some norm of the difference function 



d{s) 



--Pn{s) 
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on [so, R] is minimized. It is usually an Lp norm: 
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\d{s)\Pw{s)ds 
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where w{s) is some weight function. Obviously, different choices of p and 
w{s) will lead to different fitting coefficients a„(A^), thus different values for 
Sj^ and Afit. Since this ambiguity has not been explored in the literature on 
ACD, we have done a detailed analysis that will be presented elsewhere [§]. 
The fitting routine dependence proved to be mild. 

If the radius of the contour R is taken to be sufficiently large, the function 
jF(s) can be approximated on |s| = i? by a large momentum expansion: 



^^ h (s) 



(2.3) 
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This expression is obtained by analytically continuing the operator product 
expansion (OPE) of J-'{s) [^, which is available for both QCD-like and walk- 
ing technicolor theories, from the deep Euclidean region s -C 0. Therefore, 
we can write 

Sn = Sn M + ^tr. 



where 
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ds 
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and approximate Sn with Sn,m- The neglected term Atr is called the trun- 
cation error. It is actually not clear whether the truncation error will go to 
zero in the limit M — >■ cx) since Eq. ( |2.3| ) may be an asymptotic series and 
not a convergent one. If it is asymptotic, taking M to be larger than an 
optimum value may worsen the approximation. However, M must be taken 
to be larger than N to suppress Atr with inverse powers of R. If the series 
is convergent, the convergence may not be particularly fast requiring a very 
large M and/or R to suppress Atr- Some of these points will be discussed in 
the following, but a more complete treatment will be given in Ref . |^ . 

If the s-dependence of the expansion coefficients in Eq. ( p.3|) is "negligibly 
weak" then the approximation can be taken one step further and the s- 
dependence of the /im(s)'s dropped, i.e. 

hm{s) ^ hm{~R) = hm. 

(One must keep in mind that this is a very dangerous approximation to 
make since the analytic structure of the integrand will be completely altered 
no matter how "weak" the s-dependence may be: the imaginary part of J^{s) 
will be reduced to derivatives of 5-functions at the origin.) Therefore, 
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where 



'S'acd = 
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and we can approximate Sn,m with Sacd- Unfortunately, this is a tautolog- 
ical statement: The only way one can conclude that the s-dependence was 
"negligibly weak" is to calculate Aac and find that it is indeed suppressed. 
However, it is possible to give the following handwaving argument that Aac 
would be generically small: The difference hm{s) — h^ in the integrand of 
Aac can be expected to be smallest in the deep Euclidean region and most 



pronounced in the vicinity of the positive real s axis. But since that is where 
the difference 1/s—pn{s) is approximately zero, the contribution of this large 
difference to Aac with be highly suppressed. Thus: 

S ~ 5'acd- 

In this approximation, the integral for S'acd "will only pick up the residues of 
the single poles inside the integration contour and we find, 

mm{N,M-l} 

5'acd = 4:71 Yl an{N)hn+i. 

ra=0 

We will call Aac the analytical continuation error. 
To summarize, the ACD technique uses the relation 

S = 5'acd + Aac + Atr + Agt, 

and assumes that all three types of error can be neglected and approximates 
S with 5'acd- Therefore, the question of whether the ACD technique is 
reliable or not can be translated to the question of whether all these errors 
are under control or not. 

In actual applications of the ACD technique, there is an additional source 
of error due to the fact that the expansion coefficients hm are only known 
approximately. This is because the /i^'s are products of Wilson coefficients 
and condensates of operators. The Wilson coefficients can only be calculated 
perturbatively, and only a very limited number of them are known. The 
condensates are even less known and only a handful of them have been de- 
termined phenomenologically. As a result, the numerical uncertainty on hm 
increases with m so increasing A^ and M does not necessarily improve the 
final estimate of 5'acd since it may be swamped by the numerical uncertainty 
of the higher order terms one must include. A common strategy in the liter- 
ature is to try to choose optimum values of A^ and M such that the fit and 
truncation errors are suppressed, but the numerical uncertainty is not yet 
too large. The procedure essentially consists of the minimization of a 'total' 
error. However, it is very difficult to define the total error since it involves 
a mathematically well defined error and, on the other hand, an error which 
refiects our ignorance. The problem itself is ill-posed and difficult to treat 
with any level of rigor. 
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Table 1: S'acd and the fit, truncation, and analytical continuation errors for 
the perturbative model function. The cutoffs are [sq, -R] = [4m^,25m^], and 
the fit routine was the least square fit. The exact value of S is l/Gvr = 0.0531. 



A^ 


M 


•S'acd 


Sn,M — SacD + ^AC 


Afit 


At. 


3 


2 


0.2930 


0.0580 


-0.0002 


-0.0048 




3 


0.2883 


0.0530 




0.0002 




4 


0.2884 


0.0532 




-0.0000 


4 


2 


0.4330 


0.0632 


-0.0001 


-0.0101 




3 


0.4203 


0.0521 




-0.0010 




4 


0.4211 


0.0532 




-0.0000 




5 


0.4211 


0.0531 




0.0000 


5 


3 


0.5731 


0.0506 


-0.0000 


0.0025 




4 


0.5759 


0.0533 




-0.0002 




5 


0.5757 


0.0531 




0.0000 




6 


0.5757 


0.0531 




-0.0000 



We will therefore ignore this last type of error in this letter and confine 
ourselves to the question: if the expansion coefficients of Eq. ( p.3|) were 
known exactly, would the ACD technique give the correct result? 

3 The Perturbative Spectral Function 

Now we investigate the reliability of the ACD technique by evaluating S'acd 
and the three errors Afit, Atr, and Aac in some toy models. First we take the 
one-loop perturbative contribution of a doublet of heavy fermions to S. The 
contribution to the vacuum polarization coming from this fermion doublet 
can readily be computed Q, and therefore the S parameter can easily be 
evaluated; its value is l/67r. 

Next, we use the ACD technique to do the same computation. The per- 
turbative spectral function is 



1 m^ /■! 
47r2 s Jo 



^pert(s) = 7Z?— /_ c^xlog 



s 
1 - x{l - x)^ 



It is analytic in the entire complex s plane except for a branch cut along the 
positive real s axis which starts from s = Am?. The imaginary part of this 
function along the cut is given by 



1 777 / A-in 

Im^pe.t(s) = ms - 4m2), P = \l • 

Air s \ s 

The first few terms of the large s-expansion of J-'pcrt{s) are given by 
-47r^J'pcrt(s) 

where x = Aw? / s. Using these expressions, we calculated Sacd, '^ac; ^tr, 
and Afif The result of our calculations for several values of A^ and M are 
shown in Table ^ The fit interval was [sq, -R] = [Am'^,25ir?], and the fit 
routine was the least square fit. 

As is evident from Table |1], the fit and truncation errors are under excel- 
lent control and Sn,m reproduces the exact value of 5* accurately already at 
N = M = 3. However, the analytic continuation error is not. For the A^ = 5 
case, for instance, Sacb is larger than the exact value by more than an order 
of magnitude. In fact, we find that S'acd and Aac diverge as A^ — ;> oo. 

We conclude that neglecting the s-dependence of the /im(s)'s fails miser- 
ably as an approximation. The reason for this can be traced to the fact that 
even though the difference 1/s — Pn{s) converges to zero within its radius 
of convergence, outside it diverges. Therefore, the handwaving argument we 
gave in the previous section was wrong: the error induced by the neglect of 
the s-dependence of the /im(s)'s may be suppressed near the real s axis, but 
it is actually enhanced away from it. 

4 The Breit-Wigner Model 

In light of our result for the perturbative spectral function, the natural ques- 
tion to ask next is whether including the s-dependence of the expansion 
coefficients, i.e. not neglecting Aac, would give an accurate estimate of 5" 
for all cases. 
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Figure 2: The Spectral Functions for the Breit-Wigner Model. The ratio 
C, = Ty/my = Fa/wia is varied from C, = 1.0 to ^ = 0.1 at 0.1 intervals 
with smaller values of ^ corresponding to more prominent peaks. The arrows 
indicate the positions of the techni-p and techni-ai peaks. 



To answer this question, we consider the following model: 

1 



^v(s) 

:Fa(s) 



•2 

^ s — m\ + i^Tx 



I 



s — my + i^/sTy ' 

1 



These functions are analytic in the complex s plane except for a branch cut 
along the positive real s axis starting from s = 0. The imaginary part along 
the cut is of the Breit-Wigner type: 



\m.J^y{s) 
ImJ^A(s) 



-/^ 



-fl 



v^F 



V 



s-mlY + sVl^^'^ 



yiFA 



[s - mi)2 + sVl 



\s . 



We impose the Weinberg sum rules [H on these functions which leads to 



the constraints^ 



/v ~ /a 

r2 2 r2 2 



JTC-i 



/yJ- V 



/iri 



where /tc is the technipion decay constant. 
The value of S for this model is given by 



S = -A-K [Jv(0) - Ta 



47r 



/? 



m 



V 



/I 



ml 



47r(l + r 



f2 



m 



V 



where r = m^jnip^ 
from QCD: 



^ We fix r and the ratio j^QJvn^ using large-A^ rescaling 



m 



V 



m 



V 



m2 
m^ 



COOSA^dA^tc, 



0.4 



ml 



m 
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so that 



S ^ 0.09iVDA^TC- 



Note that the value of S does not depend on the value of the width to 
mass ratio 

ry ^ Fa 

my ITIa' 

However, the spectral function and the s dependence of the large momentum 
expansion coefficients do. In Fig. |^, we show the ^-dependence of ImjF(s). 
In the limit ^ -^ +0, these functions reduce to the well known (5-function 
(vector-meson dominance) model, and the large momentum expansion coef- 
ficients have no s-dependence. In the other limit where ^ ~ 1, we obtain a 
smooth function reminiscent of the perturbative case, and the s-dependence 
of the expansion coefficients is large. In the following, we fix ^ to the large- A^ 
value 

^ = ly ~ liL ~ 0.2, 
my irip 
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Table 2: The S'-parameter for the Breit-Wigner model function. The cutoffs 
are Sq = 0.2my and R = 5my. The dimensionless parameters are r = 0.4, 
^ = 0.2. The exact value of S is 0.257. (Since the S'ir = 0.016 correction is 
not included here, the results should be compared with 0.241). 



N 


M 


'S'acd 


Sn,m = Sacb + Aac 


Afit 


Atr 


3 


2 


0.657 


0.557 


-0.057 


-0.244 




3 


0.365 


0.102 




0.211 




4 


0.375 


0.343 




-0.030 




5 


0.375 


0.305 




0.009 


4 


2 


0.838 


0.836 


-0.006 


-0.573 




3 


-0.063 


-0.445 




0.708 




4 


0.009 


0.624 




-0.361 




5 


0.118 


0.183 




0.080 




6 


0.118 


0.244 




0.019 


5 


3 


-1.083 


-1.578 


0.014 


1.821 




4 


-0.811 


1.604 




-1.361 




5 


0.093 


0.402 




0.645 




6 


-0.230 


0.327 




-0.085 




7 


-0.230 


0.238 




0.005 


8 


7 


11.687 


-6.149 


0.000 


6.406 




8 


-3.111 


2.397 




-2.141 




10 


-1.195 


0.028 




0.234 



which is in between the two extremes. 

The large momentum expansion of our model function is given by 



:f{s) = ^ E 

^ n=0 



Itc ^ (-i)"m2r 



X2„(r)f/2„(e/2) - ^^X2„+i(r)f/2„+i(e/2) 



^Actually, the integral for the second Weinberg sum rule does not converge for this 
model. We therefore impose the equivalent constraint that the 1/Q^ term vanishes in the 
OPEfor J^(-Q2). 
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where X„ 



is defined as 



Xnir) 



while Un{cos9) = sin(n + l)^/sin^ are the Chebyshev polynomials of the 
second kind. Note that in the limit ^ = 0, the s-dependence of the coefficients 
vanishes because f/2„+i(0) = 0. 

We have a slight problem with the choice of fit interval since the branch 
cut of this model begins at s = but the IR cutoff sq must be kept non- 
zero. We will therefore choose sq to be well below the peak of the techni-p 
resonance but non-zero. The approximate form of the spectral function below 
So is 

U-^cX-3{r)V~s 



-lmjF(s) 



4-41ev^ 



171' 



TTIK 



This can be used to calculate analytically the contribution from the interval 

[0,So]: 

-S^iRlSo) = -4 / — ImJ^(s) ^ 3 >L — . 

JO s my 

The ACD technique is employed for estimating the contribution coming from 
the interval [sq, R]. We must modify the large momentum expansion of J-'{s) 
slightly to subtract out the contribution of the interval [0, Sq]: 



:F{s) =^ ^(s) - (5^ir(s) 



where 



SJ'mi-Q' 



1 r^o Im:F(s) 

n Jo s + Q^ 



E 

m=l 



[m+\)Q 



2m '■ 



and our estimate of S will be given by 

S ^ 5*^(50) + 'S'acd( So, -R)- 

For our present analysis, we used the values so/rriy = 0.2, and R/m\ = 5. 

Using these numbers and expressions, we calculated S'acd and the three 
errors Aac, Atr, and Agt for the case A^d = 1, Arc = 3. The results are 
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shown on Table 0. In this case, not only did the S'acd approximation fail, 
but even the inclusion of the s-dependence did not give accurate and stable 
results. The reason is that the truncation error is oscillatory and converges 
too slowly. It will eventually converge since the large momentum expansion 
is convergent for this model, but not before M becomes impractically large; 
in a realistic application one would have to truncate the OPE at M = 3 or 
4. 



5 Conclusions 

We have shown that the ACD technique employed by Sundrum and Hsu 
in estimating the S'-parameter is not completely reliable. We find that the 
neglect of the s-dependence of OPE coefficients is not a good approximation. 
In an example that we have analyzed we find that the truncation error is also 
not controllable, even if the spectral function is exactly known. This brings 
into doubt the utility of the ACD method in estimating the S parameter in 
complicated theories like walking technicolor. 
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